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ABSTRACT 

We simulate Shor's algorithm on an Ising spin quantum computer. The influence of 
non- resonant effects is analyzed in detail. It is shown that our "27rfc" -method suc- 
cessfully suppresses non-resonant effects even for relatively large values of the Rabi 
frequency. 
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I. The Quantum Shor Algorithm 



The quantum Shor algorithm provides an exciting opportunity for prime-factorization 
of large integers - a problem beyond the capabilities of today's powerful digital com- 
puters. Shor's algorithm utilizes two quantum registers (the x- and y- registers), which 
contain quantum bits two-level quantum systems called qubits llp-0. First, the quan- 
tum computer produces the uniform superposition of all states in the x-register - all 
possible values of x. Second, the quantum computer computes the periodic function: 
y{x) = q^imod N), where N is the number to factorize, and q is any number which is 
coprime to A^. Third, the quantum computer creates a discrete Fourier transform of 
the x-register. The measurement of the state of the x-register yields the period, T, of 
the function y{x), which is used to produce a factor of the number A^. 

In Dirac notation, the wave function of the quantum computer can be represented 
as a superposition of digital states, 

|aL_iaL_2...aiao, bM~ibM-2---bibo), (1) 

where (0 < < L — 1) denotes the state of the k-th qubit in the x-register, and 6„ 
(0 < n < M — 1) denotes the state of the n-th qubit in the ^/-register. For example, if 
the k-th qubit of the x-register is in the ground state, then: = 0, and if it is in the 
excited state, = 1. 

In decimal notation, the digital state can be represented as \x,y), where 

L-l M-l 

x=J2 «fe2^ y=Y. bnT. (2) 

fc=0 n=0 

In this notation, the initial wave function of a quantum system is: |0, 0). The uniform 
superposition of the states created in the x-register can be written as, 

*=^i:k.o). (3) 



where D = 2^ is the number of states in the x-register. After computation of the 
function y{x), we have, 

^ = ^El^'?/(^))- (4) 

V T. 



After the discrete Fourier transform, one measures the state of the x-register. The 
probability of the measurement, P(x), must be a peaked distribution with peak sepa- 
ration. Ax, equal to a multiple of 1/T. In particular, if the number of states, in the 
x-register is divisible by the period, T, then: Ax = D/T. From the value Ax one can 
find the period, T. A factor of the number can be found by computing the greatest 
common divisor of (g^^^ -f 1) and N ^ or (g^''^ — 1) and (for even T). 
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It was shown in ||^, that the simplest demonstration of the quantum Shor's algo- 
rithm can be done with only four qubits. (Two qubits represent the x-register and two 
qubits represent the ?/-register.) This primitive quantum computer is able to find a 
factor of the number 4. For = 4, the only coprime number is g = 3. The function, 
y{x) = 3^ {mod 4), has only two values: y = 1 and y = 3, and a period is, T = 2. In 
Dirac notation, the wave function (4) has the form, 

^ = i(|00,01) + |01,ll) + |10,01) + |ll,ll)). (5) 

After the Discrete Fourier transform, the wave function of the quantum computer is, 

^ = ^(|00,01) + |00,11) + |10,01) + |10,11)). (6) 

Measuring x can produce two values: x = and x = 2. So, Ax = 2, and the period is: 
P = D/ Ax = 2. Finally, g^/^ — 1 = 3 — 1 = 2, and the factor of 4 can be found as the 
greatest common divisor of 4 and 2. 

2. The Ising Spin Quantum Computer 

The Ising spin quantum computer, first introduced in ^ , consists of a one- dimensional 
chain of 1/2 spins, connected by Ising interactions. The quantum logic gates can 
be implemented in the Ising spin quantum computer using selective electro-magnetic 
pulses which induce transitions between the ground state and the excited state of a 
chosen spin. The Hamiltonian of the quantum computer during the action of the n-th 
electro-magnetic pulse can be written as, 

^ = - E [^kPk + l^kn (^e-(-"*+^")4' + e^(-"*+^")/+^ |- (7) 

s-2 
k=0 

where Uk is the resonant frequency of the k-th spin; Qkn is the Rabi frequency of 
the k-th spin during the n-th pulse; tu" and ipn are the frequency and the phase of 
the n-th electro-magnetic pulse; Jk,k+i is the constant of the Ising interaction, J| and 
= ± il^ are the operators for the k-th 1/2 spin; S is the total number of spins 
(qubits) in the quantum computer, and we choose units in which h = 1. The difference 
of resonant frequencies between spins (qubits) can be provided by differences in the 
gyromagnetic ratio of the spins or by a non-uniform external magnetic field. The 
Ising interactions also influence the resonant frequency, providing an opportunity for 
conditional logic gates. The main disadvantage of the Ising spin quantum computer 
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is the presence of non-resonant effects: a selective electro-magnetic pulse which is 
resonant to a specific spin, influences all other spins of the quantum computer. 

Destructive non-resonant effects can be significantly weakened by using a 2'n'k- 
method suggested in (see also 0], Chapter 22). The main idea of this method is 
that parameters of the electro-magnetic pulse should be chosen in such a way that a 
non-resonant spin rotates about the effective magnetic field by the angle of 2TTk (where 
k is an integer). So, at the end of the pulse, it returns to its initial position. 

To the best of our knowledge, this paper provides the first simulation of the quan- 
tum Shor's algorithm, taking into consideration non-resonant effects. We investigate 
the destructive influence of non-resonant effects and show that one can significantly 
reduce non-resonant effects by using 27rA;-pulses. 

3. Simulation of Quantum Computation 

We describe the 4-qubit quantum computer with the Hamiltonian (7) putting 3 = 4, 
uJk+i = ^^k + At^j ^kn = ^, Jk,k+i = J ■ We count spins (qubits) from the right to the 
left. The values of parameters will be given below. To implement the Shor's algorithm 
we consider application of 16 selected resonant pulses. The Hamiltonian (7) suggests 
that the n-th pulse is "cut" from a continuous harmonic oscillation, exp(za;"'t) and its 
phase is shifted by ipn- For example, the first pulse has frequency, = uj2 + 2 J, and 
phase, (pi = n/2; the second pulse has frequency, u"^ = + J, and phase, = tt/2, 
and so on. We will use single- integer decimal notation for the 4-spin basic states, i.e., 

\aiao, bibo) = \p), p = bo + 2bi + 2^ao + 2^ai, 

and Cp is the amplitude corresponding to the state \p). The Schrodinger equation for 
the amplitude, Cp, of the state, can be written as, 

15 

= c^Vp^e^^[i{Ep- E^)t + rp^{uj'^t + ipn)\. (8) 

Here we use the interaction representation, i.e. Cp — > Cpexp{—iEpt); Vpm is the matrix 
element between the states \p) and \m) which is equal to Vpm = —fl/2 for single- 
spin transitions and zero for other transitions; Ep and E^, are the energies of the 
corresponding states, r^^ = —1 for Ep > E^,, and Vpm = 1 for Ep < Em- If one 
neglects non-resonant effects, the ra-th pulse produces only a rotation of a resonant 
spin which can be described as a transformation of the basic states: 

|0) ^cos(a„/2)|0) + «exp(-Z¥P„)sin(«„/2)|l), (9) 

|1) cos(a„/2)|l) + iexp(iv9„) sin(a„/2)|0), 

where a„ = fir^ is the angle of rotation of the average spin, and r„ is the duration of 
the n-th pulse. For example, first two pulses are the 7r/2-pulses, i.e. ai = 0^2 = vr/2. 
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Instead of a direct discrete Fourier transform we have used an idea of Coppersmith 
and Deutsch, first described in ||^. Following this idea, the final wave function (without 
non-resonant effects) contains four states: 

|00,01), |00,11), |01,01), |01,11). (10) 

One must reverse the results of the measurement of the x-register to get the actual 
result of the discrete Fourier transform. Following this rule, one obtains: x = and 
X = 2 with equal probability, 1/2. 

It follows from (10) that in the resonant approximation we have, 

|Ci|2=|c3p=|c5p=|c7|2 = l/4, (11) 

and zero probabilities for all other states. 

4. Results of Simulations of the Shor's Algorithm 



Fig. la, shows the actual values of the "resonant" probabilities (11) for the following 
values of parameters, 

Acj = 10, J = 1, n = 0.1. (12) 

One can see that the expression (11) is approximately satisfied. Small deviations of 
the resonant probabilities from 1/4 are associated with nonzero probabilities for "non- 
resonant" states. Fig. lb, (with a magnified scale) shows that non-resonant states are 
very non-uniformly excited: the probability of the upper state, |15), is close to 10~^; 
the probabilities of the states |9) and |11) are close to 10~^; and all other non-resonant 
probabilities have smaller values. Fig. 2 shows a radical change of the results for a small 
increase of the Rabi frequency (from f2 = 0.1 to f2 = 0.112). Now the probabilities 
of the resonant states are significantly different: the probability of measuring x = 2 is 
noticeably greater than for the value x = 0. The probabilities of error connected with 
non-resonant states has sharply increased. 

Fig. 3 shows the return to the quasi-resonant picture for Q = 0.125: the proba- 
bilities of the resonant states (9) are close to 1/4, and the probabilities of the non- 
resonant states are less than 10~^. Figs 4 through 7 demonstrate the same periodicity 
at increasing the Rabi frequency: "quasi-ideal" implementation of Shor's algorithm 
at Q = 0.1666, Q = 0.25 (Figs 5 and), and its destruction for intermediate values, 
Q = 0.1458, Q = 0.2083 (Figs 4 and 6). Fig. 8 demonstrates the dependence of reso- 
nant probabilities on the value of the Rabi frequency, fl, which continuously approaches 
the value of the Ising constant, J. The resonant probabilities are close to (11) at fixed 
values of fl: 

n = 0.1, 0.125, 0.1666, 0.25. (13) 
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The last point where the probabihties approach the value 1/4 is: n — 0.51639, but the 
deviation from 1/4 at this point is much greater than for previous values of Q (13). 
Fig. 9 demonstrates the destruction of the Shor's algorithm when the value of the 
frequency difference between neighboring qubits, Aa;, approaches the value of the Ising 
constant, J. (The Rabi frequency, Q and J are given by (12).) Significant deviation of 
the probabilities from their "resonant" values begins at approximately Au = 3. The 
greatest probability of error occurs in the upper state, |15). Fig. 10 shows the time 
dependence of the resonant probability, |c3p, during the action of the last two pulses, 
for different values of Aa;. One can see, that for Aa; > 4, the function |c3(t)p roughly 
does not depend on Aa;. For Aa; < 3, this function depends significantly on Aa;. 

Note that, in Figs 1-10, the angles of rotation, did not change, i.e. the increase 
of fl was compensated by the decrease of the pulse dutation, r„. Fig. 11 shows the 
destruction of Shor's algorithm when the phases, <^„, deviate randomly from their cor- 
rect values, <^°. The phases were chosen randomly in the interval (</?° — e,ip'^ + s), 
for four values of e. Significant influence of the phase fluctuations can be observed at 
£ > 0.5 rad. Again, the main error occurs in the upper state, |15). Fig. 12 shows the 
influence of the fluctuation of the angle of rotation, a„ on the performance of Shor's 
algorithm. We have randomly changed the values of r„ in the interval (r° — e,T^ + e), 
for different values of e. The values of fl, Aa;, and r° are given by (12), where T^is the 
correct value of r^. The noticeable destruction in the performance of Shor's algorithm 
was observed at £ > 2. The corresponding angle of rotation is: Qe > 0.2 rad. 

5. Discussion 

The most important result of our simulation of Shor's algorithm is the successful pe- 
riodic suppression of non-resonant effects when the value of Rabi frequency, fl, ap- 
proaches the Ising constant, J. We shall now discuss this phenomenon in detail. 

Suppose the n-th electro-magnetic pulse induces a non-resonant transition between 
the states \p) and |m), where Ep > Em- The equation of motion (8) for these two 
states can be written approximately as, 

icp = -{n/2)cmQw{i[{Ep -Em- uj^^Y - (14) 

icm = -{ft/2)cpQ-x.Y){i[{Em - Ep + u'^)t + (/?„]}. 

Here we assume that the probability of the non-resonant transition between the states 
\p) and \m) is much greater than any other transition probability involving these two 
states. With an accuracy up to the phase, the solution of (14) can be written as, 

Cp{tn) = Cp(t„_i)[cOs(fieTn/2) - {iA/Qe) Sm{Qern/'2)] + C„i(t„_i) [(iil/Qe) sin(neT"n/2)] , 

(14a) 

Cm{tn) = C„(tn-l)[cOs(QeT„/2) + {iA/ft^) sin(QeTn/2)] + Cp{tn-l)[{ift/fte) Sin(QeT„/2)]. 
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Here is the effective frequency for the non-resonant transition, 



(15) 



and A = Ep — — uj"', (t„_i, t„) is the time-interval of the n-th pulse; t„ — = x„. 

One can see that for QeTn = 27r/c [k = 1,2,...) the probability of non-resonant 
transition vanishes. This is the main idea of the "27r/c" method 0, For a vr-pulse 
{flTn = vr), the values of Q which satisfy this 27r/;;-condition are, 

Q = |A|/V4A;2 - 1 = |A|/V3, \A\/Vl5, |A|/V35, |A|/V63, ... (16) 

For a 7r/2-pulse, the corresponding values of Q are, 

Q = |A|/Vl6A;2 - 1 = \A\/Vl5, \A\/V63, ... (17) 

If the Rabi frequency satisfies the 27r/c-condition (17) for a 7r/2-pulse, it automatically 
satisfies the condition (16) for a vr-pulse. 

Let us consider the most probable non-resonant transitions in our system. We 
assume that J -C Au, so the most probable non-resonant transitions occur when the 
frequency difference, |A|, is of the order of J. For the right spin {k = 0) and the left 
spin {k = 3) we have two frequencies depending on the state of the only neighbor: 
ujk + J and Uk — J- If the r/ pulse is turned on either of these transitions, the frequency 
difference, |A|, for the other one will be 2J. For the inner spins [k = 1,2), we have 
three possible resonant frequencies: Uk, Uk — 2J and uJk + 2J. The frequency difference, 
|A|, can take two values: 2 J and 4 J. 

One cannot suppress both non-resonant transitions with |A| = 2 J and |A| = 4 J. 
Putting |A| = 2J in (17) we obtain, 

n = 2J/V16P - 1 ^ J/2k. 

For |A| = 4 J, we have, 

n = 4J/Vl6A;2 - 1 ^ J/k. 

So, if one satisfies the 27rA;-condition for |A| = 2 J, one satisfies approximately the 
same condition for |A| = 4 J. But total suppression of two non- resonant transitions is 
impossible. 

Now we consider an example - the suppression of non-resonant effects during the 
first two pulses. The initial state of our system is the ground state, |00, 00). First, we 
apply a 7r/2 pulse with the frequency, uj2 + 2J, which is resonant with the spin with 
the number k = 2. The non- resonant effects for three other spins are small if Au is 
large enough. After the action of the first pulse, we have the superposition, 

(l/y2)(|00,00) + |01,00)). (18) 
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Next, we apply a 7r/2-pulse with the frequency uj-^ + J. This pulse is resonant with the 
= 3 spin (in its ground state, the first term in (18)). But this pulse also perturbs the 
same spin in the excited state (the second term in (18)). The frequency of the transition 
for the spin with = 3 in the second term in (18) is uj^ — J. So, the frequency difference 
is: |A| = 2J. Taking a value of Vt which satisfies (17) with |A| = 2J, one suppresses 
this non-resonant transition. 

Putting f2 = 2 J/ VlS/c^ — 1 with J = 1 we obtain the values: 

VL ^ 0.100, 0.125, 0.167, 0.252, 0.516, 

which are close to the values which were obtained from the computer simulations. 

To decrease the time of quantum computation, one should use the largest possible 
values of fl. If the acceptable probability of error state is of the order of 10"'^, then 
one should use the value Q ~ 0.25 (Fig. 7). If a probability of the error state of 10~^ is 
acceptable, one can use the value, Q ~ 0.516. The time of quantum computation can 
be further reduced if one uses different values of Q for 7r/2- and vr- pulses. According 
to (16) and (17), the maximum value of the Rabi frequency, Q, for a 27r rotation of 
non-resonant spins is, Q ~ 0.51639 for 7r/2-pulses, and fl ~ 1.1547 for vr-pulses. Fig. 
13 shows the probabilities of the states for this minimal-time Shor's algorithm imple- 
mentation. The greatest probability of error ~ 10~^ is connected with the states |11) 
and 1 15). Note, that for 7r-pulses the value of Rabi frequency is greater than J and the 
27rA;-method still suppresses the non-resonant effects. 

Conclusion 

We reported the first simulation of the quantum Shor's algorithm for prime factor- 
ization with taking into consideration non-resonant effects. We have considered the 
Ising spin quantum computer of four qubits: two qubits in the x-register, and two 
qubits in the ^/-register. This primitive quantum computer is able to factor the small- 
est composite number, = 4. While a simulation of much more complicated system 
of 15 qubits was reported in , that work deals with well-separated ions in an ion trap 
and it does not take into consideration non-resonant effects. 

We have studied the destructive influence of non-resonant effects when the Rabi 
frequency, Q, approaches the value of the Ising interaction constant, J. We also studied 
the influence of the random fluctuations of the duration and the phase of the electro- 
magnetic pulses. 

The main results of our consideration are the following: 

1. When the Rabi frequency, Q, approaches the Ising interaction constant, J, there 
is no monotonic increase of destructive non-resonant effects. Non-resonant effects are 
suppressed effectively when the value of Q satisfies, 

n = 2J/Vl6A;2 -1, k= 1,2,.... (19) 
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for all pulses. These values correspond to 27rk rotation of a non-resonant spin whose 
frequency is the most close to the resonant frequency. This method of suppression of 
non- resonant effects (27r/c-method) was first suggested in |^ and discussed in P, ^, ||. 
But this method has never been verified for a quantum computer algorithm. 

2. A 27r/i;-method allows one to decrease the total time of computation by taking 
the maximum possible value of Q. The minimum time of computation was achieved at 
the values: f2 = 2J/\/lE for 7r/2-pulses, and = 2J/\/3 for vr-pulses. The probability 
to measure an error state at these values of Q was of the order 10~^ (Fig- 13). 

3. The probability of error caused by non-resonant effects or by random fluctuation 
of parameters is connected with the relatively small amount of states. As a rule, one, 
two, or three error states have the probability of excitation much greater than all 
other error states taken together. Normally, the upper state has a large probability of 
excitation. 

The next step should include bigger number of qubits, and application of 2Tck- 
method together with error correction approaches. 
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Fig. 1: Probabilities of the states |n), |c„p, (n = 0, . . . , 15) at the end of the Shor's 
algorithm. The parameters are: Ao; = 10, J = 1.0 and Q = 0.1; (a) the probabilities 
for the "resonant" states, with complex amplitudes, ci, C3, C5 and cy; (b) the probabil- 
ities for other states (which remain small), at different scale. 



Fig. 


2: 


The same as in Fig. 1, for Q, = 0.112. 


Fig. 


3: 


The same as in Fig. 1, for Q — 0.125. 


Fig. 


4: 


The same as in Fig. 1, for Q = 0.1458. 


Fig. 


5: 


The same as in Fig. 1, for fl = 0.1666. 


Fig. 


6: 


The same as in Fig. 1, for Q = 0.2083. 


Fig. 


7: 


The same as in Fig. 1, for fl = 0.2083. 


Fig. 


8: 


Probabilities, ci p, cap, C5 p and C7 p at the end of the Shor's algorithm, as 



a function of the Rabi frequency, Q, . The parameters are: Acu — 10 and J = 1.0. 

Fig. 9: Probabilities, |c„p, (n = 0, . . . , 15) at the end of the Shor's algorithm, for 
different values of Acu. The parameters are: fl = 0.1 and J = 1.0. 

Fig. 10: Dynamics of the probability, |c3(i)|^, near the end of the Shor's algorithm, 
for different values of Auj. The parameters are: J — 1 and Q, — 0.1. Arrows show the 
end of 14-th, 15-th, and 16-th pulses. 

Fig. 11: Probabilities, |cnp, (n = 0, . . . , 15) at the end of the Shor's algorithm, for ran- 
dom phases: (pk = -\- Aip^'^ <^^°^ are the required phases of the pulses {k = 1, 16); 
A(pk are random phases which were chosen in the interval: A(p e [— e, e]; (1) e = 0, (2) 
e = 0.1, (3) e = 0.5, (4) e = 0.8. The parameters are: Auj = 10, Q = 0.1 and J = 1.0. 

Fig. 12: Probabilities, |cnp, (n = 0, . . . , 15) at the end of the Shor's algorithm, 
for random duration of resonant pulses, r^, (n = 1,...,16); = r^^ + At„, where 
is a required duration of the n-th pulse, and At„ is a random variable uniformly 
distributed between — e and +€. The parameters are: Alo = 10, Q = 0.1 and J — 1.0. 

Fig. 13: Probabilities, |cnp, n = 0, . . . , 15, at the end of the Shor algorithm, for 
four qubits. The parameters are: Au; — 10, J — 1.0, Q — 2/-\/3 for 7r-pulses, and 
n = 2/VlE for 7r/2-pulses. 
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